Search for: All records

Excited-state properties of crystalline organic semiconductors are key to organic electronic device applications. Machine learning (ML) models capable of predicting these properties could significantly accelerate materials discovery. We use the sure-independence-screening-and-sparsifying-operator (SISSO) ML algorithm to generate models to predict the first singlet excitation energy, which corresponds to the optical gap, the first triplet excitation energy, the singlet–triplet gap, and the singlet exciton binding energy of organic molecular crystals. To train the models we use the “PAH101” dataset of many-body perturbation theory calculations within the GW approximation and Bethe–Salpeter equation (GW+BSE) for 101 crystals of polycyclic aromatic hydrocarbons (PAHs). The best performing SISSO models yield predictions within about 0.2 eV of the GW+BSE reference values. SISSO models are selected based on considerations of accuracy and computational cost to construct materials screening workflows for each property. The screening targets are chosen to demonstrate typical use-cases relevant for organic electronic devices. We show that the workflows based on SISSO models can effectively screen out most of the materials that are not of interest and significantly reduce the number of candidates selected for further evaluation using computationally expensive excited-state theory. 

Abstract The excited-state properties of molecular crystals are important for applications in organic electronic devices. TheGWapproximation and Bethe-Salpeter equation (GW+BSE) is the state-of-the-art method for calculating the excited-state properties of crystalline solids with periodic boundary conditions. We present the PAH101 dataset ofGW+BSE calculations for 101 molecular crystals of polycyclic aromatic hydrocarbons (PAHs) with up to  ~500 atoms in the unit cell. To the best of our knowledge, this is the firstGW+BSE dataset for molecular crystals. The data records include theGWquasiparticle band structure, the fundamental band gap, the static dielectric constant, the first singlet exciton energy (optical gap), the first triplet exciton energy, the dielectric function, and optical absorption spectra for light polarized along the three lattice vectors. The dataset can be used to (i) discover materials with desired electronic/optical properties, (ii) identify correlations between DFT andGW+BSE quantities, and (iii) train machine learned models to help in materials discovery efforts. 

The excited-state properties of molecular crystals are important for applications in organic electronic devices. The GW approximation and Bethe-Salpeter equation (GW+BSE) is the state-of-the-art method for calculating the excited-state properties of crystalline solids with periodic boundary conditions. We present the PAH101 dataset of GW +BSE calculations for 101 molecular crystals of polycyclic aromatic hydrocarbons (PAHs) with up to ∼500 atoms in the unit cell. The data records include the GW quasiparticle band structure, the fundamental band gap, the static dielectric constant, the first singlet exciton energy (optical gap), the first triplet exciton energy, the dielectric function, and optical absorption spectra for light polarized along the three lattice vectors. In addition, the dataset includes the density functional theory (DFT) single-molecule and crystal features used in Liu et al. [npj Computational Materials, 8, 70 (2022)]. We envision the dataset being used to (i) identify correlations between DFT and GW +BSE quantities, (ii) discover materials with desired electronic/ optical properties in the dataset itself, and (iii) train machine-learned models to help in materials discovery efforts. We provide examples to illustrate these three use cases. 

The excited-state properties of molecular crystals are important for applications in organic electronic devices. The GW approximation and Bethe-Salpeter equation (GW +BSE) is the state-of-the-art method for calculating the excited-state properties of crystalline solids with periodic boundary conditions. We present the PAH101 dataset of GW +BSE calculations for 101 molecular crystals of polycyclic aromatic hydrocarbons (PAHs) with up to ∼500 atoms in the unit cell. The data records include the GW quasiparticle band structure, the fundamental band gap, the static dielectric constant, the first singlet exciton energy (optical gap), the first triplet exciton energy, the dielectric function, and optical absorption spectra for light polarized along the three lattice vectors. In addition, the dataset includes the density functional theory (DFT) single-molecule and crystal features used in Liu et al. [npj Computational Materials, 8, 70 (2022)]. We envision the dataset being used to (i) identify correlations between DFT and GW +BSE quantities, (ii) discover materials with desired electronic/ optical properties in the dataset itself, and (iii) train machine-learned models to help in materials discovery efforts. We provide examples to illustrate these three use cases. 

Integrating second order nonlinear (χ(2)) optical materials on chip is an ongoing challenge for Si photonics. Noncentrosymmetric molecular crystals have the potential to deliver high χ(2) nonlinearity with good thermal stability, but so far have been limited to growth from solution or the melt, which are both difficult to control and scale up in manufacturing. Here, we show that large (>100 μm) single crystal domains of the nonlinear molecule 2-[3-(4-hydroxystyryl)-5,5-dimethylcyclohex-2-enylidene] malononitrile (OH1) can be grown monolithically on either glass or Si via vacuum evaporation, followed by a short thermal annealing step. The crystallites are tens of nanometer thick and exhibit strong second harmonic generation with their primary χ(2) tensor component lying predominantly in plane. Remarkably, we find that a single domain can grow uninterrupted through nearby channels etched on a Si wafer, which may provide a path to integrate OH1 on Si or Si3N4 waveguides for a broad range of χ(2)-based photonic integrated circuit functionality. 

Thermally activated delayed fluorescence (TADF) is the internal conversion of triplet excitons into singlet excitons via reverse intersystem crossing (RISC). It improves the efficiency of OLEDs by enabling the harvesting of nonradiative triplet excitons. Multiple resonance (MR) induced TADF chromophores exhibit an additional advantage of high color purity due to their rigid conformation. However, owing to the strict design rules there is a limited number of known MR-TADF chromophores. For applications in full-color high-resolution OLED displays, it is desirable to extend the variety of available chromophores and their color range. We computationally explore the effect of chemical modification on the properties of the MR-TADF chromophore quinolino[3,2,1-de]acridine-5,9-dione (QAD). QAD derivatives are evaluated based on several metrics: The formation energy is associated with the ease of synthesis; The spatial distribution of the frontier orbitals indicates whether a compound remains an MR-TADF chromophore or turns into a donor-acceptor TADF chromophore; The change of the singlet excitation energy compared to the parent compound corresponds to the change in color; The energy difference between the lowest singlet and triplet states corresponds to the barrier to RISC; The reorganization energy is associated with the color purity. Based on these metrics, QAD-6CN is predicted to be a promising MR-TADF chromophore with a cyan hue. This demonstrates that computer simulations may aid the design of new MR-TADF chromophores by chemical modification. 

The true molecular conformation and the crystal structure of benzo[ e ]dinaphtho[2,3- a ;1′,2′,3′,4′- ghi ]fluoranthene, 7,14-diphenylnaphtho[1,2,3,4- cde ]bisanthene and 7,16-diphenylnaphtho[1,2,3,4- cde ]helianthrene were determined ab initio by 3D electron diffraction. All three molecules are remarkable polycyclic aromatic hydrocarbons. The molecular conformation of two of these compounds could not be determined via classical spectroscopic methods due to the large size of the molecule and the occurrence of multiple and reciprocally connected aromatic rings. The molecular structure of the third molecule was previously considered provisional. These compounds were isolated as by-products in the synthesis of similar products and were at the same time nanocrystalline and available only in very limited amounts. 3D electron diffraction data, taken from submicrometric single crystals, allowed for direct ab initio structure solution and the unbiased determination of the internal molecular conformation. Detailed synthetic routes and spectroscopic analyses are also discussed. Based on many-body perturbation theory simulations, benzo[ e ]dinaphtho[2,3- a ;1′,2′,3′,4′- ghi ]fluoranthene may be a promising candidate for triplet–triplet annihilation and 7,14-diphenylnaphtho[1,2,3,4- cde ]bisanthene may be a promising candidate for intermolecular singlet fission in the solid state. 

Abstract Singlet fission (SF), the conversion of one singlet exciton into two triplet excitons, could significantly enhance solar cell efficiency. Molecular crystals that undergo SF are scarce. Computational exploration may accelerate the discovery of SF materials. However, many-body perturbation theory (MBPT) calculations of the excitonic properties of molecular crystals are impractical for large-scale materials screening. We use the sure-independence-screening-and-sparsifying-operator (SISSO) machine-learning algorithm to generate computationally efficient models that can predict the MBPT thermodynamic driving force for SF for a dataset of 101 polycyclic aromatic hydrocarbons (PAH101). SISSO generates models by iteratively combining physical primary features. The best models are selected by linear regression with cross-validation. The SISSO models successfully predict the SF driving force with errors below 0.2 eV. Based on the cost, accuracy, and classification performance of SISSO models, we propose a hierarchical materials screening workflow. Three potential SF candidates are found in the PAH101 set.